Capillary  Discharge  Thruster 
Experiments  and  Modeling 


Robert  Martin1 


ERC  Inc.1,  In-Space  Propulsion  Branch, 
Air  Force  Research  Laboratory 
Edwards  Air  Force  Base,  CA  USA 


Electric  propulsion  systems 
June  2016,  Rhode-Saint-Genese,  Belgium 


Distribution  A:  Approved  for  public  release;  Distribution  Unlimited;  PA#  16279 


& 


^  'O  O 


U.&  AIR 


<  □ 


Outline 


O  Introduction 
Q  AFRL  CDT  Experiments 
Q  CDT  and  Related  Models 
o  Current  Status  &  Future  Work 
e  Conclusion 


:  □  I  g  >  <  ^  ►  <  -f  ;  JF  'OQlO' 


Distribution  A:  Approved  for  public  release;  Distribution  Unlimited;  PA 


R.S.  Martin  (ERC  Inc.) 


16279 


2/30 


Electric  Propulsion 


Spacecraft  Propulsion  Relevant  Plasma: 

o  From  hall  thrusters  to  plumes  and 
fluxes  on  components 

©  Complex  reaction  physics  i.e. 
Discharge  and  Breakdown  in  FRC 

©  Relevant  Densities  often  Span 
6+  Orders  of  Magnitude 

©  Spatial  scales  of  interest  span 
jim-lOOm  range 
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Electric  Propulsion  Plumes 


FRC 
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Chamber  Environment 


Electric  Propulsion 


Spacecraft  Propulsion  Relevant  Plasma: 

o  From  hall  thrusters  to  plumes  and 
fluxes  on  components 

o  Complex  reaction  physics  i.e. 
Discharge  and  Breakdown  in  FRC 

o  Relevant  Densities  often  Span 
6+  Orders  of  Magnitude 

©  Spatial  scales  of  interest  span 
/im-l00m  range 

All  Relatively  Low  Energy  Density... 
Connection  to  HEDP  Capillary  Discharge? 


All  Rarefied  Gas  and  Plasma 


...and  Highly  Non-Equilibrium 


& 


<  □ 


'O  Qs  O' 


Interest  in  Pulsed 


Need  Efficiency  Across  Thrust  Range: 
o  Spacecraft  Power  is  Constrained 
o  Fundamental  Tradeoff:  Isp  vs.  Thrust 

o  Optimal  Trade  Mission  Dependent 
(i.e.  Station  Keeping  vs.  Orbit  Insertion) 


Thrusters 
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Interest  in  Pulsed 


Need  Efficiency  Across  Thrust  Range: 
o  Spacecraft  Power  is  Constrained 
o  Fundamental  Tradeoff:  Isp  vs.  Thrust 

o  Optimal  Trade  Mission  Dependent 
(i.e.  Station  Keeping  vs.  Orbit  Insertion) 

o  Electrothermal  -  Electrostatic  Gap 


Thrusters 


SPECIFIC  IMPULSE,  SEC 

Burton,  et.  al„  AIAA  Paper,  (TDS83-10),  1983. 
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Interest  in  Pulsed 


Need  Efficiency  Across  Thrust  Range: 
o  Spacecraft  Power  is  Constrained 
o  Fundamental  Tradeoff:  Isp  vs.  Thrust 

o  Optimal  Trade  Mission  Dependent 
(i.e.  Station  Keeping  vs.  Orbit  Insertion) 

o  Electrothermal  -  Electrostatic  Gap 

o  Burton  Predicted1 /Demonstrated2 
Efficient  Pulsed  Electrothermal  (PET) 

burton,  Goldstein,  Tidman,  Winsor,  AIAA  Paper,  82(1952),  1982. 

2Burton,  Fleischer,  Goldstein,  Tidman,  Winsor,  NASA,  (CR-179464),  1984. 
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Motivation 


Capillary  Discharge  Thruster  Viability: 

o  Efficiency  Gap  for  Moderate  ISP  EP 
(750s-3000s) 

o  Capillary  Discharge  Conceptually 
Efficient  (rjt> 65%)  in  this  Range 

o  Burton  Predicted  r]t  ^80% 

o  Burton  Observed  only  56%  Max 
(0.085  Ns  @  1600s  Isp) 

«  Realizing  full  Efficiency  requires 
Deeper  HED  Physics  Understanding 

e  CDTs  are  Simple  Small  Devices 
Accessible  to  Lab  Experiments 

Must  Converge... 

Theory,  Models,  and  Experiments 
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Burton,  et.  al.,  AIAA  Paper,  (TDS83-10),  1983. 


Capillary  Discharge 


Key  Processes  for  Design  &  Efficiency: 

©  Assumptions 

-Unmagnetized/LTE/Coupling 

©  Energy  Balance 

-Conduct/Evaporate/Bond  Break/Ionize 

©  Ablation 

-Photo-ablation/Macro-particles/Pyrolysis 

©  Radiative  Transport 

-Optical  Depth/Spectrum 

©  Ionization/Recombination 


-Breakdown/Recombination  Rate 


'ROCESS 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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CD  Process: 


Key  Processes  for  Design  &  Efficiency: 

©  Magnetization 

-Essentially  Electrothermal 
-Weaker  Assumption  if  n=(9(le24/m3) 


Plasma-/?: 


R  =  Ez 
P  Pb 

Pj  —  nkT  &  Pb 


D  _  M q/ 

27 rr 


£2 

2/^o 


£ 


8tz&T  /7rr\2 
MO  W/ 


Using: 

r= 2eV,  n=le25/m3,  r=2mm,  I=6kA 


/3=22»1 


CD  Process:  Assumptioi 


Key  Processes  for  Design  &  Efficiency: 
©  Magnetization 

©  Local  Thermodynamic  Equilibrium 
-LTE:  Acceleration/Collision  Balance 
-Highly  Collisional  After  Breakdown 


LTE  Parameter: 

K=^eSr 
A  <*<->,  =  r(2j) 

if  _  1  c6  me  (  n 

A  ”  128  7r2e^3  m-  UW 

Using: 

T= 2eV,  n=le25/m3,  E- 1  e5  V/m 


CD  Process:  Assumptioi 


Key  Processes  for  Design  &  Efficiency: 
o  Magnetization 

o  Local  Thermodynamic  Equilibrium 

o  Plasma  Coupling 

-Potential/Kinetic  Energy  Balance 
-Degree  Ideal  Plasma  EOS  Applies 
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Non-Ideal  Parameter: 

P  _  UPE 

Uke 


UpE  47reo  r 

UKE  =  kT 


e2n1/3 

4neo 


e2n1/3 

47T6()kT 


r  = 

Using: 

T= 2eV,  n=le25/m3 


T^O.16 

Non-Ideal?  T  <  1,  but  not 

r«  1? 

Ideal  Assumption  Used,  but 
Should  be  Revisited. 
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Key  Processes  for  Design  &  Efficiency: 


o  Energy  Balance 

-Evaporate  C2H4  fron  wall:  0.5eV 

-Break  C-C  Bond:  4.5ev 

-Break  4  C-H  Bond: _ 14.0eV 

Total  6  Atoms:  19.0ev 

-Dissociation  Energy/Atom:  3.2eV 

-Mean  Ionization  Energy/ Atom:  12.8ev 

Total  Energy/Ion:  16.0ev 


Would  be  Energy  Sink  Inhibiting  Efficiency...  but 
Recombination  before  Exit  Captures  Ion  Energy! 


Pancotti,  PhD  Dissertation,  USC,  2009. 


Losses  via  Radiation/Conduction  to  Housing... 

Limited  on  Discharge  Timescales 
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CD  Process:  Ablation 


Key  Processes  for  Design  &  Efficiency: 

©  Photo-Ablation 

-Direct  Ablation  by  Photon  Energy 
-Polymers  Highly  Susceptible  to 
Photo-Ablation 

-Still  Requires  <300nm  Photons 
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CD  Process:  Ablation 


Key  Processes  for  Design  &  Efficiency: 

©  Photo-Ablation 

-Direct  Ablation  by  Photon  Energy 
-Polymers  Highly  Susceptible  to 
Photo-Ablation 

-Still  Requires  <300nm  Photons 

©  Macro-Particles 

-Subsurface  Energy  Deposition 
-Vaporaztion  Ejects  Macro-particles 
-Particles  Evaporate  Quickly  (S/V  ratio) 
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CD  Process:  Ablation 


Key  Processes  for  Design  &  Efficiency: 

©  Photo-Ablation 

-Direct  Ablation  by  Photon  Energy 
-Polymers  Highly  Susceptible  to 
Photo-Ablation 

-Still  Requires  <300nm  Photons 

©  Macro-Particles 

-Subsurface  Energy  Deposition 
-Vaporaztion  Ejects  Macro-particles 
-Particles  Evaporate  Quickly  (S/V  ratio) 

©  Pyrolysis 

-Thermal  Evaporation 

-Surface  Temperature  must  Exceed  Tvap 

-Conductive  Losses  with  Tvap 
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Key  Processes  for  Design  &  Efficiency: 


o  Photo-Ablation 

-Direct  Ablation  by  Photon  Energy 
-Polymers  Highly  Susceptible  to 
Photo-Ablation 

-Still  Requires  <300nm  Photons 

e  Macro-Particles 

-Subsurface  Energy  Deposition 
-Vaporaztion  Ejects  Macro-particles 
-Particles  Evaporate  Quickly  (S/V  ratio) 

o  Pyrolysis 

-Thermal  Evaporation 

-Surface  Temperature  must  Exceed  Tvap 

-Conductive  Losses  with  Tvap 


Direct  Ablation  Preferable...  Spectrum? 


CD  Process:  Radiation 


Key  Processes  for  Design  &  Efficiency: 
o  Spectrum 

-Several  Materials  Investigated 
-Spectra  Generated  using  PrismSpect® 

o  Optical  Depth 

-A$»0(1)*-0(O.1)J? 

-High  Radiation  Conductivity  — >>  Uniform  T 


100  200  300  400  500  600  700  800  900  1000 

Wavelegth  (nm) 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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CD  Process:Ionization 


Key  Processes  for  Design  &  Efficiency: 


o  Ionization  Process 
-Wire  Breakdown  is  Chaotic 


CD  Process:Ionization 


Key  Processes  for  Design  &  Efficiency: 

o  Ionization  Process 
-Wire  Breakdown  is  Chaotic 
-Paschen  Breakdown  more  Predictable 
-Breakdown  Voltage  Density  Dependent 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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CD  Process:Ionization 


Key  Processes  for  Design  &  Efficiency: 

o  Ionization  Process 
-Wire  Breakdown  is  Chaotic 
-Paschen  Breakdown  more  Predictable 
-Breakdown  Voltage  Density  Dependent 
-Spark 
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CD  Process:Ionization 


Key  Processes  for  Design  &  Efficiency: 

®  Ionization  Process 
-Wire  Breakdown  is  Chaotic 
-Paschen  Breakdown  more  Predictable 
-Breakdown  Voltage  Density  Dependent 
-Spark 

o  Recombination 
-Recombination  Rate: 
ve  =  a3n2  =  8.75  x  10“277’-9/2n2Hz 
(T  in  eV,  n  in  cm- 3) 

-Thermal  Velocity  /  Mean  Free  Path: 


Device 

Burton  PET 

Density,  n* 
Temp,  T* 
Rate,  ve 
Velocity,  u* 
MFP,  A* 

5.4e27/nr5 

4ev 

5.0el4  Hz 

1 .2e4  m/s 
2.4e-ll  m 

Device 

Pancotti  CDT(Est) 

Density,  n* 
Temp,  T* 
Rate,  ve 
Velocity,  u* 
MFP,  A* 

1 .0e25/mJ 

2ev 

3.9el0  Hz 

9.6e3  m/s 

2.5e-7  m 

2.5e-7  m<c2mm 
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Experimental  Setup 


Thruster  Design  &  Ignition: 


©  Baseline 


Housing 


Polyethylene  Tube 


1 —  — 

|- - UU(W - 

- m' - j 

l  J 

X  Charging  Circuit  — 

1 — — 

_  Discharge  Circuit 

- m - 

\  Dump  Circuit  ^ 

£  Capillary  plasma  forms  a  resistive  load 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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Experimental  Setup 


Thruster  Design  &  Ignition: 
o  Baseline 

o  Wire  Ignition 

-  Simple  and  Reliable 

-  Chaotic  Process 

-  Random  Ignition  Delays 

-  Bi-Modal  Performance 

-  Only  Single  Use 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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Experimental  Setup 


Thruster  Design  &  Ignition: 
o  Baseline 

o  Wire  Ignition 

-  Simple  and  Reliable 

-  Chaotic  Process 

-  Random  Ignition  Delays 

-  Bi-Modal  Performance 

-  Only  Single  Use 

©  Paschen  Ignition 

-  More  Repeatable 

-  Enabled  Use  of  Thrust  Stand 

-  Requires  some  Background  Density 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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Experimental  Setup 


Thruster  Design  &  Ignition: 
o  Baseline 

o  Wire  Ignition 

-  Simple  and  Reliable 

-  Chaotic  Process 

-  Random  Ignition  Delays 

-  Bi-Modal  Performance 

-  Only  Single  Use 

©  Paschen  Ignition 

-  More  Repeatable 

-  Enabled  Use  of  Thrust  Stand 

-  Requires  some  Background  Density 

o  3 -Electrode  Ignition 

-  More  Applicable  to  Space  Vacuum 

-  Dielectric  Flashover 

-  Less  Electrode  Erosion 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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Experimental  Facility 


& 


□ 


Torsional  Thrust  Stand 


70  (t)  +  CQ  ( t )  +  K@  ( t )  =  M  ( t ) 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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Torsional  Thrust  Stand 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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Paschen  Ignition 


Key  Processes  for  Design  &  Efficiency: 


o  Linear  Impulse  with  Energy 
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Paschen  Ignition 


Key  Processes  for  Design  &  Efficiency: 
o  Linear  Impulse  with  Energy 
o  Large  Scatter  in  Isp 
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Paschen  Ignition 


Key  Processes  for  Design  &  Efficiency: 
0  Linear  Impulse  with  Energy 
0  Large  Scatter  in  Isp 
0  Large  Scatter  in  Efficiency 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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Paschen  Ignition  Performance 


Key  Processes  for  Design  &  Efficiency: 

0  Linear  Impulse  with  Energy 

0  Large  Scatter  in  Isp 

0  Large  Scatter  in  Efficiency 

o  Performace: 

-Impulse:  20-100  mNs 
-Specific  Impulse:  350-700s 
-Efficiency:  8-17%  (Nozzleless  Design) 


Scatter  due  to  Electrode  Erosion? 


Energy  (J) 

Pancotti,  PhD  Dissertation,  USC,  2009. 
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Paschen  Ignition 


Key  Processes  for  Design  &  Efficiency: 
o  Same  Linear  Impulse  with  Energy 
o  Better  Correlated 


RMANCE 


0.12 


Energy  (J) 
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Paschen  Ignition 


Key  Processes  for  Design  &  Efficiency: 
©  Same  Linear  Impulse  with  Energy 
o  Better  Correlated 
©  Clearer  Isp  Trends 
©  Higher  Isp  when  Longer 
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Paschen  Ignition  Performance 


Key  Processes  for  Design  &  Efficiency: 
©  Same  Linear  Impulse  with  Energy 
©  Better  Correlated 
©  Clearer  Isp  Trends 
©  Higher  Isp  when  Longer 
©  Longer  also  Higher  Efficiency 
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Paschen  Ignition  Performance 


Key  Processes  for  Design  &  Efficiency: 

©  Same  Linear  Impulse  with  Energy 

o  Better  Correlated 

o  Clearer  Isp  Trends 

©  Higher  Isp  when  Longer 

©  Longer  also  Higher  Efficiency 

©  8cm  Efficiency  Constant  with  Energy 

©  Performace: 

-Impulse:  30-120  mNs 
-Specific  Impulse:  350-650s 
-Efficiency:  9-17%  (Nozzleless  Design) 


0.12 


Energy  (J) 


Pancotti  PhD  Dissertation  USC  2009 


Temperature  from 


Spitzer  Resistivity: 

o  Ratio  of  Rate  Electron  Momentum 
Exchange  to  Current  Density: 
r]=l/a  =  ^f 
o  For  a  Lorentz  Gas: 


9 


(Stationary  Ions/Nonineracting  Electrons) 

_ pi3l2Zm2ec2vlnA 

^ L  ~  2{2kT)3/2 

With  e-e  Collsions  (Spitzer-Harm) 


Ionic  Charge  Z 

1 

2 

4 

16 

oo 

7 E 

0.582 

0.683 

0.785 

0.923 

1.000 

Distribution  A:  Approved  for  pi 
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Optical  Diagnostics 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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Optical  Emission  Spectrc 


Time  Resolved  OES: 

o  Uses  Spectral  Line  Shape  not  Intensity 

-Simpler  Calibration 

o  Pulsed  Device  Requires  Time  Resolved 

o  Kinetics  Mode  via  Pixel  Time  Shifts 
-5pixel/Spectra 
-16 /a  s/Spectra 
-O.lnm/Pixel 
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Optical  Emission  Spectrometry 


Time  Resolved  OES: 

o  Uses  Spectral  Line  Shape  not  Intensity 

-Simpler  Calibration 

o  Pulsed  Device  Requires  Time  Resolved 


Spectra  5 

r 

□ 

□ 

□ 

□ 

Spectra  4 

£ 

s 

s 

s 

£ 

Spectra  3 

s 

s 

s 

£ 

Spectra  2 

Is 

gg 

Spectra  1 

o  Kinetics  Mode  via  Pixel  Time  Shifts 
-5pixel/Spectra 
-16/x  s/Spectra 
-O.lnm/Pixel 
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Electron  Density  Diagnostic 


Hydrogen-a  Line  Broadening: 

9  Neutral  Broadening 
(Small) 

9  Doppler  Broadening 

(Ha,  H^:<lA) 

9  Resonance  Broadening 
(N-N  Collisions) 

9  Van  der  Waals  Broadening 
(Also  N-N  Collisions) 

9  Stark  Broadening 
O(10nm) 


Af/2  «  1  x  10“4[A] 

Af/2  =  7.16xlO-7A0^[A] 

Af/2  =  8.6  x  lO-30^fkX20XrfrNi[A\ 

Af/2  »  3  x  KT^cf  (J)3/I0N[A] 
Aj\2  «  2.5  x  10-9a1/2N2/3[A] 
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Electron 


Hydrogen-a  Line  Broadening: 

9  Neutral  Broadening 

9  Doppler  Broadening 

9  Resonance  Broadening 

9  Van  der  Waals  Broadening 

9  Stark  Broadening 


Distribution 
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Diagnostic 


Huddlestone  &  Leonard,  Plasma  Diagnostic  Techniques,  Academic  Press,  ’65. 


~  2.5  x  10-* * * * 9a1/2Ne2/3[A] 
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Electron  Density  Diagnostic 


Hydrogen-a  Line  Broadening: 

9  Neutral  Broadening 

9  Doppler  Broadening 

9  Resonance  Broadening 

9  Van  der  Waals  Broadening 

9  Stark  Broadening 


Huddlestone  &  Leonard,  Plasma  Diagnostic  Techniques,  Academic  Press,  ’65. 
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& 


:  □ 


'O  Qv  O' 


Optical  Emission  Spectrometry 


Time  Resolved  Electron  Density: 
o  Spectrum  fit  to  Lorentzian  Profile: 

/(A  -  A0)  =  £  [^]  where  27=FWHM 


t  2 


Pancotti, 


Optical  Emission 


Time  Resolved  Electron  Density: 
o  Spectrum  fit  to  Lorentzian  Profile: 

/(A  -  Ao)  =  £  Where  27=FWHM 

o  Fit  inverted  for  vs.  Time 
-Density  le25/m3  Estimate 
-Optical  Depth  — w  Exit  Plane  nel 


^OMETRY 


ne  vs.  Time,  6cm  Capillary 


ne  vs.  Time,  2500V  Discharge 

Pancotti,  PhD  Dissertation,  USC,  2009. 
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OP  Slab  Model: 


Conservation  of  Mass: 

2AW  ■  r  —  AeneCi 


T/  _  dn  _ 

V  ’  dt  — 


-  AeneCevh  -  2AW@ 


o  Conservation  of  Energy: 

t/  d(ne)  y  I1  / 

V  '  ~dT  V  '  ~o\n^T) 

Where: 

n  is  the  Plasma/Gas  Density 

V  is  the  Slab  Volume 

Aw  is  the  Wall  Area 

Ae  is  the  Exit  Area 

T  is  the  Ablation  Flux 

Cs  is  Sound  Speed  (at  the  Exit) 

Superscript- ()e  is  Isentropically  Expanded  Exit  Value 


e  is  the  Energy  Density 
I  is  Current 

a(n,  T)  is  the  Conductivity 
h  is  the  Enthalpy 
Cs  is  the  Sound  Speed 
0  is  the  Wall  Energy  Flux 


Pekker,  40th  AIAA  Plasmadynamics  and  Laser  Conference,  2009. 
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Modeling:  OD 


OD  Slab  Model: 


a  Conservation  of  Mass: 

V  ■  f  =  2AW  •  T  -  AeneCes 

o  Conservation  of  Energy: 

V-^  =  V-  ^  -  AeneCesh  -  2 Aw@ 

Slab  Capillary,  La  =  8  cm,  Da  =  4  mm  Slab  Capillary,  La  =  8  cm,  Da  =  4  mm 


Pekker,  40th  AIAA  Plasmadynamics  and  Laser  Conference,  2009. 
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Modeling:  ID 


1DPDE  Model: 

o  Conservation  of  Mass: 

^ r  +  m  [(apu )1  =  A* 


+  m  Ka(pu+p ))]  =P^~A™f 


Conservation  of  Momentum: 

d(Apu)  ,  0 

dt 

Conservation  of  Energy: 

¥  +  §~x  \{Au{e  +p))}  =  A  (i Qj  -  Qrad  -  Qconv  -  Qab  -  <f>) 


Where: 

A  is  the  Cross  Section  Area 
p  is  the  Mass  Density 
u  is  the  Velocity 
Aw  is  the  Wall  Surface  Area 
r  is  the  Ablation  Mass  Flux 

p  is  the  Pressure 
/  is  the  Viscous  Wall  Friction 

Pancotti,  PhD  Dissertation,  USC,  2009. 


e  is  the  Total  Energy 

Qj  is  the  Joule  Heating 

Qrad  is  the  Radiant  Energy  Losses 

Qconv  is  the  Convection  Energy  Losses 

Qab  is  the  Ablation  Energy 

<F  is  Viscous  Wall  Energy  Loss 
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9  Conservation  of  Mass: 

+  m  \tAPu)\  =  Aw  ■ r 

o  Conservation  of  Momentum: 

^Itr^  +  miiMpu+p))]  =P^~Awf 

o  Conservation  of  Energy: 

^r~  +  §^c  \{Au{e  +  p))]  =  A  ( Qj  —  Qrad  —  Qconv  —  Qab  — 

Energy  Flux  for  5cm  2500V  Capillary  Discharge  Comparison  of  ID  Model-ne  and  Experimental  Observation 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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Comparison  of 


Discharge  Current  Predictions: 

o  Comparison  of  2500V 
Discharge 

©  Similar  Profiles/Trends 

©  Wire  Highest  Current 

©  Paschen  Lowest  Current 

©  Models  Over-Predict  Tail 
(Especially  0D) 

©  Initial  dl/dt  Incorrect 
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and  Experiment 


Pancotti,  PhD  Dissertation,  USC,  2009. 
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-3 -Electrode  System  developed  by  Pancotti 
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-Nozzle  added  for  Efficient  Energy  Conversion 
-Additional  Propellant  Materials  were  Studied 
-High  Efficiencies  Demonstrated,  Isp~  1000s 
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Additional  Thruster  Development 


For  Proof-of-Concept  Demonstration: 

©  Repeatable  Ignition 

-3 -Electrode  System  developed  by  Pancotti 

o  Desired  Isp  &  77 

-Nozzle  added  for  Efficient  Energy  Conversion 
-Additional  Propellant  Materials  were  Studied 
-High  Efficiencies  Demonstrated,  Isp~  1000s 

©  Robust  Propellant  feed  Mechanism 
-Remains  Unresolved 
-Burton  studied  Liquid/Gas  Schemes 
-Additional  research  Required 


& 


<  □ 


'O  A  O' 


Extension  to  Particle  Kinetic  Models 


Breakdown  non-LTE: 

o  Many  Particles  —^Continuous  Distribution 


& 


-<  □ 


'O  c*s  O' 


Extension  to  Particle  Kinetic  Models 


<  □ 


'O  c*s  O' 


Extension  to  Particle  Kinetic  Models 


Breakdown  non-LTE: 

o  Many  Particles  —^Continuous  Distribution 

o  Discretized  VDF  yields  Vlasov  Models 
But  3D3V  Severe  Dimensionality  Curse 

o  Particle  Methods  (i.e.  PIC)  Simplify  to  Delta 
Functions 

o  Sparse  Representation  in  3D3V  Phase  Space 
But  Added  Noise  and  Low  Dynamic  Range 


VDF 


& 


<  □ 


Extension  to  Particle 


Breakdown  non-LTE: 

o  Many  Particles  —^Continuous  Distribution 

o  Discretized  VDF  yields  Vlasov  Models 
But  3D3V  Severe  Dimensionality  Curse 

o  Particle  Methods  (i.e.  PIC)  Simplify  to  Delta 
Functions 

o  Sparse  Representation  in  3D3V  Phase  Space 
But  Added  Noise  and  Low  Dynamic  Range 

®  Collisions  between  Discrete  Velocities 
o  Inelastic  Collisions  in  Tail  Impact  High  Moments 


etic  Models 


VDF 

(Av)2  >  l„ 

I  I  III  □ 


& 


<  □ 


'O  c*s  O' 


Extension  to  Particle  Kinetic  Models 


Breakdown  non-LTE: 

o  Many  Particles  —^Continuous  Distribution 

o  Discretized  VDF  yields  Vlasov  Models 
But  3D3V  Severe  Dimensionality  Curse 

o  Particle  Methods  (i.e.  PIC)  Simplify  to  Delta 
Functions 

o  Sparse  Representation  in  3D3V  Phase  Space 
But  Added  Noise  and  Low  Dynamic  Range 

o  Collisions  between  Discrete  Velocities 
o  Inelastic  Collisions  in  Tail  Impact  High  Moments 
®  All-or-Nothing  Collision  — >  Rare  Large  Events 


VDF 

(Av)2  >  lD 


-/ 


4 


& 


-<  □ 


'O  c*s  O' 


Extension  to  Particle 


Breakdown  non-LTE: 

o  Many  Particles  —^Continuous  Distribution 

o  Discretized  VDF  yields  Vlasov  Models 
But  3D3V  Severe  Dimensionality  Curse 

o  Particle  Methods  (i.e.  PIC)  Simplify  to  Delta 
Functions 

o  Sparse  Representation  in  3D3V  Phase  Space 
But  Added  Noise  and  Low  Dynamic  Range 

®  Collisions  between  Discrete  Velocities 
o  Inelastic  Collisions  in  Tail  Impact  High  Moments 
o  All-or-Nothing  Collision  — >  Rare  Large  Events 
o  Variable  Weights  — >•  more  DOF  representing  Tails 


Extension  to  Particle  Kinetic  Models 


Breakdown  non-LTE: 

o  Many  Particles  —^Continuous  Distribution 

o  Discretized  VDF  yields  Vlasov  Models 
But  3D3V  Severe  Dimensionality  Curse 

o  Particle  Methods  (i.e.  PIC)  Simplify  to  Delta 
Functions 

o  Sparse  Representation  in  3D3V  Phase  Space 
But  Added  Noise  and  Low  Dynamic  Range 

®  Collisions  between  Discrete  Velocities 
o  Inelastic  Collisions  in  Tail  Impact  High  Moments 
o  All-or-Nothing  Collision  — >  Rare  Large  Events 
©  Variable  Weights  — >•  more  DOF  representing  Tails 
©  Fractional  Collisions  New  Numerical  Particles 


4 


Extension  to  Particle  Kinetic  Models 


Breakdown  non-LTE: 

o  Many  Particles  —^Continuous  Distribution 

o  Discretized  VDF  yields  Vlasov  Models 
But  3D3V  Severe  Dimensionality  Curse 

o  Particle  Methods  (i.e.  PIC)  Simplify  to  Delta 
Functions 

o  Sparse  Representation  in  3D3V  Phase  Space 
But  Added  Noise  and  Low  Dynamic  Range 

®  Collisions  between  Discrete  Velocities 
o  Inelastic  Collisions  in  Tail  Impact  High  Moments 
o  All-or-Nothing  Collision  — >  Rare  Large  Events 
©  Variable  Weights  — >•  more  DOF  representing  Tails 
©  Fractional  Collisions  New  Numerical  Particles 
©  Conservative  Merge  Needed  to  Control  Growth 


VDF 


□  •  <  g 


Distribution  A:  Approved  for  public  release;  Distribution  Unlimited;  PA 


'O  c*s  O' 


R.S.  Martin  (ERC  Inc.) 


*l#  16279 


26/30 


Extension  to  Particle  Kinetic  Models 


Breakdown  non-LTE: 

o  Many  Particles  —^Continuous  Distribution 

o  Discretized  VDF  yields  Vlasov  Models 
But  3D3V  Severe  Dimensionality  Curse 

o  Particle  Methods  (i.e.  PIC)  Simplify  to  Delta 
Functions 

o  Sparse  Representation  in  3D3V  Phase  Space 
But  Added  Noise  and  Low  Dynamic  Range 

®  Collisions  between  Discrete  Velocities 
o  Inelastic  Collisions  in  Tail  Impact  High  Moments 
o  All-or-Nothing  Collision  — >  Rare  Large  Events 
©  Variable  Weights  — >•  more  DOF  representing  Tails 
©  Fractional  Collisions  New  Numerical  Particles 
©  Conservative  Merge  Needed  to  Control  Growth 

4  □  4  gP  ►  4  B  ►  <5  :  ^ 


Distribution  A:  Approved  for  public  release;  Distribution  Unlimited;  PA 


R.S.  Martin  (ERC  Inc.) 


L#  16279 


26/30 


Octree  Binning 


Phase-Space  Decomposition 


©  Given  a  Set  of  Particles... 
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Phase-Space  Decomposition 
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©  Particles  Binned  in  Octants 
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Octree  Binning 


Phase-Space  Decomposition 


o  Given  a  Set  of  Particles... 

o  Particles  Binned  in  Octants 

o  Octants  Recursively  Sub-Divided 

o  Recursion  Halted  at  1 -Particle/Bin  or 
Other  Criteria  such  as  Bin-Density 

Restricts  Phase-Space  Diffusion  to 
Within  Local  Bins 

(Entropy,  ^  n  log(n).  Constant  within  Octree  Adaptive  Quadrature  ) 
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ID IV  Particle  Kinetic  Breakdown  Model 


250V  DC-Diode  Test  Case: 
o  Full  3D  Electrostatic-PIC 
o  Averaged  to  ID  XT-Plot 
o  250V  Cathode  — )>  Anode 
Q  MCC-Ionization  Collisions 
o  Secondary  Emission  at  Cathode 


Number  Density,  (#/m3)  Number  Density,  (#/m3) 
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ID IV  Particle  Kinetic  Model-Breakdown 


IKY  DC-Diode  Test  Case: 

0  Voltage  increased  to  1KV 
0  Otherwise  Identical  to  250V 
0  Much  Stronger  Ionization 


Position,  (mm) 
Number  Density,  (#/m3) 


Control  Merged 
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